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Abstract 
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1 Introduction 

The study of the Riccati differential equations goes back to the early days of modern mathe¬ 
matical analysis, since such equations represent one of the simplest types of nonlinear ordinary 
differential equations and consequently Riccati equations play an important role in physics, 
mathematics and engineering sciences. The usual first-order Riccati equation appears as a 
reduction from a second-order linear ordinary differential equation when taking into account 
invariance under dilations according to Lie recipe (the inverse property is called the Cole-Hopf 
transformation). Such correspondence between first-order Riccati equations and second-order 
linear ODEs can be beautifully manifested through their solutions. For example, if one solu¬ 
tion of a linear second-order ODE is known, then it is known that any other solution can be 
obtained by means of a quadrature (i.e. a simple integration), and the same property holds 
true for the Riccati equation, but with two quadratures. Moreover, first-order Riccati equa¬ 
tion is a prototypical example of (systems of) differential equations admitting a superposition 
rule, also called Lie systems [ml [13]. This is why some authors [23] have considered it as the 
first step in the study of nonlinear (systems of) differential equations. 

The second-order Riccati equation also becomes of paramount importance in recent years, 
mainly because of its connection with the integrable ODEs. It is known that the second-order 
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Riccati equations are related to a large class of integrable ordinary differential equations of an- 
harmonic oscillators, viz. the Ermakov-Pinney, Normalised Ermakov-Pinney, Neumann type 
systems, etc. [H |2ll [32l |36l ESI 09]. It has been studied in [H] from a geometric perspective 
and it has been proved to admit two alternative Lagrangian formulations, both Lagrangians 
being of a non-natural class (neither potential nor kinetic term). A more geometric approach 
can be found in [S] where the theory of Darboux polynomials, the extended Prelle-Singer 
methods, (pre-)symplectic forms, or Jacobi multipliers, are used. 

Higher-order Riccati equations can be obtained by reduction from the Matrix Riccati 
equation [10l|57]. Not only hrst-order Riccati equation but also higher-order Riccati differen¬ 
tial equations can be linearised via Cole-Hopf transformation (i.e. they appear as a reduction 
from linear differential equations when taking into account the dilation symmetry). 

The study of higher-order Riccati equations [8] was also carried out from the perspective 
of the theory of Darboux polynomials and the extended Prelle-Singer methods. Higher-order 
Riccati equations play the role of Backlund transformations for integrable partial differential 
equations of order higher than that of the KdV equation. In fact Grundland and Levi [HI] 
constructed the Backlund transformations of several integrable systems from higher-order 
Riccati equations. Such Backlund transformations are closely related to many important 
integrability properties such as the inverse scattering method, the Painleve property. Lax 
pairs, and an inhnite number of conservation laws [69] . 

In the study of differential equations one hnds cases that are in some sense solvable, or 
integrable, which enable one to study their dynamical behaviours using Lie theoretic methods. 
They often admit different geometric and Hamiltonian formulations. Integrable systems are a 
fundamental class of explicitly solvable dynamical systems of current interest in mathematics 
and physics. One notable example is the Milne-Pinney equation |6ll9l[IIl|53l|6l]: in 1950 
Pinney presented in a one page paper [6T] the solution of the equation 

y" + n‘^{x)y = 4 , 

yS 

where the symbol y" denotes the second derivative with respect to the independent variable 
X. He gave the solution in the form 

y{x) = {A(pl + 2B 01 02 + O 

where yi = (j)i{x) and y 2 = 02 (a;) are any two linearly independent solutions of the associated 
linear equation y" + fl‘^{x)y = 0, and A, B and C are related according to — AC = l/W"^ 
with W being the constant Wronskian of the two linearly independent solutions. Ermakov 
introduced the above equation as an additional auxiliary equation to a second order linear 
differential equation to dehne a system of two second order differential equations, and found 
an invariant when multiplying by an integrating factor and obtained the invariant after inte¬ 
gration with respect to x. 

The solutions of Painleve’s Second Equation (or Painleve H) 

y" = a + xy + 2y^ (1) 
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are meromorphic functions in the plane |67]. It is known (see e.g. the paper by Gromak 
[29] and [30]) that every transcendental solution w has infinitely many poles with residue +1 
and also with residue —1, except when a = ±1/2 and w solves the Riccati equation (see [67] i 


dw 

dx 


±-±w 

2 
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In the paper we derive the Painleve II, as well the second-order second degree Painleve II 
[221 03] using Riccati hierarchy, which is in turn connected to the Virasoro stabilizer set or 
projective connection. The second-order second degree Painleve II is also known as Jimbo- 
Miwa equation |15] and plays an important role in random matrix theory [26] . 


1.1 Motivation and plan 

Our main goal is to study higher-order Riccati equations in various directions. In particular, 
we focus onto the geometrical aspects of the Riccati sequence. At hrst we explore its con¬ 
nection to the stabilizer set of the Virasoro orbit, also known as the projective vector field 
equation. There are many papers devoted to Virasoro algebra and projective connection on 
S^, but very few papers discuss about their relations to nonlinear oscillator equations, Riccati 
chains and (second-order) Painleve equations. This paper provides a relationship between the 
stabilizer set of the Virasoro orbit and second-order Riccati equations. Many papers (for ex¬ 
ample, |33l|33[59]) have been written to elucidate the interaction between the geodesic flows 
on the Bott-Virasoro group (see for example, [13 HHl EO]) and integrable partial differential 
equations, especially the KdV type systems, but very few articles deal with finite-dimensional 
systems. 

Some of the main questions to be discussed in this paper can be summarized in the 
following points: 

• Relation of Riccati equation with the Virasoro algebra. 

The Virasoro group plays a very important role in integrable systems. It is known that 
the Virasoro group serves as the conhguration space of the KdV and the Camassa-Holm 
equations and these equations can be regarded as equations of the geodesic flows related 
to different right-invariant metrics on this group. In other words, they have the same 
symmetry group. Recall that the Virasoro group is a one-dimensional central extension 
of the group of smooth transformations of the circle. In this work, we focus on the 
integrable dynamical systems connected to the stabilizer set of the Virasoro orbit. In 
fact the Milne-Pinney equation, which describes the time-evolution of an isotonic oscil¬ 
lator (also called pseudo-oscillator) ~ i.e. an oscillator with inverse quadratic potential- 
is the first and foremost example of this class. It is known [331 EZ] that the entire 
(coupled) KdV family is connected to the Euler-Poincare formalism of the (extended) 
Bott-Virasoro group. This connection can be extended to the super Bott-Virasoro group 

PI. 

The second-order Riccati equation has also a nice geometric interpretation in terms of 
Virasoro orbit [3S1ES]- Vector helds f{x)d/dx G X{S^) associated to the stabilizer set 
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of Virasoro orbit are called projective vector fields and the corresponding differential 
equation satisfied by the function / is called the projective vector held equation [M1I35] . 
It is known that the second-order Riccati equations are associated to these equations. In 
fact, solutions of a large class of 0 -|- 1 dimensional integrable systems can be expressed 
in terms of the global and local projective vector helds. 

• Analysis of the inhnitesimal symmetries of the standard Riccati equation. 

This study is presented by making use of a geometric approach. It is proved that these 
symmetries are related to the projective vector held equation. 

• Study of the relation of the second-order Riccati equation with Painleve equations. 

There are a certain number of properties relating the second-order Riccati equation with 
some other equations of the mathematical physics as the Airy equation or the Milne- 
Pinney equation [53l|6l]. Moreover, the second-order Riccati equation is also related to 
some of the PainlevAGambier equations and this can be used to hnd the solutions of 
Painleve’s Second Equation (or Painleve II) for special values of parameters. 

The paper is organised as follows. Section 2 is a short presentation of properties of hrst- 
and higher-order Riccati equations and their relations with linear equations and nonlinear 
superposition rules. In Section 3 we give a concise introduction to Virasoro orbit, stabilizer 
set and projective connections and in Section 4 we study the inhnitesimal symmetries of 
standard Riccati equation and their relations with the projective vector held equation. We 
elucidate in Section 5 the connection between the ordinary Painleve II and second-order 
second degree Painleve II equations with the Virasoro orbit in Section 5. 


2 First-order and higher-order Riccati equations 


Our main goal is to explore the higher-order Riccati equations, but before giving the formal 
description of higher-order Riccati equations and Riccati chain, a little introduction may be 
helpful. So we start our journey by giving a short introduction to hrst-order Riccati equation. 
The usual Riccati equation 

u'= f{x) + g{x)u + h{x)u^, (2) 


is a hrst-order nonlinear diherential equation with a quadratic non-linearity. The solutions of 
the Riccati equation are free from movable branch points and can have only movable poles 
[SI mi E]. This is the simplest nonlinear diherential equation admitting a superposition 
rule for expressing the general solution in terms of particular solutions. The linear fractional 
transformation (or change of variables) 


a{x)u -|- b{x) 
c{x)u -|- d{x) ’ 


ad — be = 1, 


(3) 


transforms each Riccati equation into another one [12]. In order to solve a Riccati equation 
by quadratures, it is enough to know one particular solution which allows us to reduce the 
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problem to carry out two quadratures, while when two particular solutions are known the 
problem can be reduced to a new differential equation solveable by just one quadrature, and 
hnally if three particular solutions ui{x), U 2 {x), u^^x), are known, we can construct all other 
solutions u without use of any further quadrature [13]. This is carried out by using the 
property that the cross-ratio of four solutions is a constant (i.e. for any other solution u{x), 
there is a real number k such that: 


u{x) — ui{x) u^{x) — ui{x) 

— fC ' 


u{x) - U2{x) U^{x) - U2{xy 


(4) 


where k is an arbitrary constant characterising each particular solution. For instance ui{x) 
corresponds to k = 0 and U 2 {x) is obtained in the limit of fc —)■ cxo). 

It was shown by Lie and Scheffers [51] that the Riccati equation is essentially the only first- 
order ordinary nonlinear differential equation of which possesses a nonlinear superposition rule 
which comes from the preceding relation: 


u = ^{ui,U2,U3; k) = 


Ul{u3 - U2) -h kU2{ux - M3) 


U3-U2 + k{ui - Ms) 

(i.e. if Ui{x),U 2 {x) and M 3 (x) are solutions of (|2D, then for any real value fc, 

Mi(a;)(M3(a;) - U2{x)) ku2{x){ui{x) - U3{x)) 


u{x) = 


Usix) - U2{x) + k{ui{x) - U3{x)) 


is another solution and each solution is of this form for an appropriate choice of k). 

Therefore the general solution u{x) is non-linearly expressed in terms of three generic 
particular solutions, Mi(a;), U2{x) and u^^x), and a constant parameter k. 

Remark also that a first-order Riccati differential equation (]2]) with h{x) of constant sign 
in an open set, transforms under the relation 


M = — 


1 ^ 

h y 


(5) 


into a second order differential equation 


y" + ai(x)y' + aQ(x)y = 0 


( 6 ) 


with 

h' 

ai = -g-—, ao = fh. (7) 

This relation ([5]) is not a change of coordinates or transformation but for any real number 
A, y and Xy have associated the same function u. Conversely, given the linear second-order 
differential equation ([6]), for any nowhere vanishing function h, (|5]) reduces (|6]) to the Riccati 
equation ([2]) with coefficients determined by (171) . 

In a similar way, consider the following nonlinear second-order equation in an open set I 

m- 


u" + [f 3 o{x) + f 3 i{x)u~\u'+ ao{x) + ai{x)u + a2{x)u^ + a3{x)u^ = 0 , x G / C M, ( 8 ) 
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where we suppose that a^^x) > 0,Wx E I and the two functions (3q, (3i, are not independent 
but satisfy 


/^o 


Oi2 



Pi = 3v^. 


In particular, the case /Sq = «! = 02 = 0, = 1, Pi = 3, was studied in Davis and 

Ince books [23l HI] and will be more carefully considered in next sections, while that of 
Po = chq = ai = a 2 = 0, Pi = 1, was studied by Leach and coworkers [l9l |50] who proved 
that only the case a^lx) = 1/9 is linearizable, possesses eight symmetries and is completely 
integrable. 

An important property is that, as indicated in [T3], the nonlinear equation ([HD can be 
transformed into a third-order linear equation 


!/"' + ai(x)y" + ax(x)y' + ai(x)y = 0 


(9) 


by the substitution 


u{x) 


1 y'{.x) 

y{x) 


( 10 ) 


The equation (|HD is a second-order generalization of the usual first-order Riccati equation 
and because of this is usually known as the Riccati equation of second-order [T3]: it is a 
nonlinear equation whose solution of can be expressed in terms of solutions of a third order 
linear differential equation. Conversely, if we start with a linear third-order equation then 
we can get the associate nonlinear second-order equation by introducing the appropriate 
reduction of order. 

One can go further and define j-order Riccati equations as those appearing as a reduction 
of a (j-l-l)-order linear differential equation [8]. More specifically, given a differential equation 
y{n) _ jnvariance under dilations suggests, according to Lie recipe, to look 

for a new variable ^ such that the dilation vector field yd/dy becomes {l/k)d/dz, the factor 
k being purely conventional (we often put k = 1). Then y = up to an irrelevant factor, 
and with this change of variable, y^^'^ = 0, for n > 1, becomes R'^~^{u) = 0 with u = z 
and R the differential operator R = {D + ku). In fact, when D = d/dx, we have that 
e-kzDe^^ = D + ku, and then we can prove by complete induction that = ke^^RJ^~^{u) 
(see [8]): the property is true for n = 1, because y' = k z = k u = ke^^ R^{u), and if 
the property is true till n = j, it also holds for n = j -t- 1, because 


D^+^y = D{y^^'>) = D{ke^^W-\u)) = ke^\D + ku)R^-\u) = ke^^W{u). 


Therefore, Lie recipe applied to the invariant under dilations equation = 0 transforms 
such an equation into R^'^~^\u) = 0. The first terms of such Riccati sequence are 


R^{u) = u, R^{u) = u'+ ku‘^, R'^{u) = u” + 3kuu'+ k'^u^, ... (11) 

R^{u) = u'"+4:uu"+3u'^+6u‘^ u'+u‘^, R'^{u) = u^^^+5uu'"+10u'u"+15uu'‘^+10u^u"+l0u^u'+u^,... 

( 12 ) 
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and for the general linear differential equation of order n, 


ao{x)y+ = 0, an{x) = 1, 

i=i 


( 13 ) 


reduces to 

n 

ao{x) + aj{x) R^~^{u) = 0 . 

i=i 

3 Virasoro algebra, projective vector field and second 
order Riccati equaltion 

Let = T*S^ be the cotangent bundle of a circle S^. Since is diffeomorphic to the 
l-dimensional Lie group U{1), there exists a canonical coordinate x, dehned up to a factor, 
which we can £x such that the domain of the chart is (0,27r). This is a trivial real line 
bundle on S^. Similarly, let denote the m-fold tensor product of The local coordinate 

expression of a section of is given by s{x) = g{x) dx'^, where dx™ is a shorthand notation 
for the tensorial product of m times dx. The set £(^2”^) of such sections is a free 
module, a natural basis being given by dx™, and in this sense such sections are described by 
functions in 5^. A vector held in is of the form 

X, = f(x)-^eX(,S^), ( 14 ) 

where / is a function on (i.e. a 27r-periodic function on the real line). Taking into account 
that 

Cxfdx^ = m f dx^ 

we see that the inhnitesimal action of Xj on a section s = g{x) dx™ of f2™ is given by (see 
e.g. [31]) 

^XfS = ifg' + mf'g)dx^, (15) 

where Cxf is the Lie derivative with respect to the vector held Xf given by (1TT|) . 

Sometimes it is convenient to use the above mentioned identihcation of a section with its 
component, a function, and then the Lie derivative when acting on sections of VP expressed 
in terms of functions, denoted R^Xfi dehned as 

= (16) 

This really means that if s = '0(x) dx^, then 

CxfS = +Pf'{x)'4){x'^ dx^ = dx^. 
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There is a natural symmetric bilinear map x —)■ given by the 

usual commutative product of tensor densities 


f{x){dxr^ ® g{x){dxr^ ^ f{x)g{x){dxr^+^\ 


In particular for m = 0, r(ff°) = C°°(S'^). We can at least formally extend the values of the 
index m to the set of integer numbers, or even to rational numbers, and then the product 
by elements of r(f2°) is the external composition law of the module structure. Note that the 
elements of when multiplied by those of r(ff™) play as an inner contraction with a 

vector held and in this sense is to be identihed with X(S'^) (i.e. we can denote the 

tangent bundle as (= TS^)). 

Note also that when Xg = g{x)d/dx G then putting p = —1 in flT6|) we get the usual 
expression: 


^XfXq - 






(17) 


We will denote the ‘square root’ of the tangent and cotangent bundle of respec¬ 

tively. The respective sections of such bundles will be of the form s = 'ip{x)dx'^^ G r(n^2). 

Another remarkable property is that the C'°°(S'^)-module r(f2™') of sections of the m-fold 
tensor product bundle fl™' has a natural structure of a Poisson algebra with the commuta¬ 
tive product given by the above usual product of tensor densities and the Poisson algebra 
commutator given by the Rankin-Cohen bracket 


{f{x){dx)^,g{x){dx)^} = {mf{x)g'{x) - nf'{x)g{x)){dx)^^"^^\ 


which is also known as hrst transvectant. 

The group Diff (S^) acts on P(ff™'), as given by 


^*{f{x){dx)"^) = fi^{x)){^'{x))^{dx)"^ 

for $ G Diff (S*^). The Lie algebra of such group is identihed with the Lie algebra X(S'^) and 
acts on this space P(D™') by the Lie derivative given in flT5|l . 

The theory of projective connections on the circle has much to do with the theory of 
Riccati equations we are dealing with. We will use the following dehnition of a projective 
connection |13]: 

Definition 1 A projective connection on the circle is a linear second-order differential oper¬ 
ator 

A : P(D-^) ^ P(Di) (18) 

such that: 

1. the principal symbol of A is the identity and 

2. [ (Asi)s 2 = [ Si(As 2), for all pairs of sections Si ET{n-'2). 

Js^ isi 
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Let us take s = 'ilj{x)dx 2 G r(f2 2 )^ then a linear second-order differential operator is 
such that As e is locally described by |13] 

As = {atp"-\-htp'-\-ctp) dx^ ^ (19) 

where a, b and c are real numbers. 

The first condition in the definition of the projective connection, implies that for A to be 
a projective connection it must be a = 1, and, on the other side, the second condition implies 
that 6 = 0. Hence each projective connection can be identified with a Hill operator 

A„ = ^ + v{x), (20) 

dx^ 

where v is an arbitrary function in (i.e a periodic function in M). 


Definition 2 Given a projective connection 


d^ 


( 21 ) 


a vector field Xf G X(S'^) is called projective vector field with respect to the projective connec¬ 
tion Ay when it leaves invariant the projective connection (i.e. 


C^'^AyS = Ay{C™s), 


’(- 1 / 2 ), 


( 22 ) 


for all s 2 )^ where is the expression for the Lie derivative with respect to Xf given 

in /fT^H. 


The characterization of projective vector fields with respect to a given projective connec¬ 
tion is given in the following theorem (see H): 

Theorem 1 Xf G r(f2“^) is a projective vector field with respect to the projective connection 
ISP if and only if f satisfies (see m) 

r + Avf + 2v'f = 0 . (23) 

-i.e. f is a solution of the third-order linear differential equation called projective vector field 
equation: 

y'" + Av y'+ 2v'y = 0. (24) 

Proof: In fact, if s = 'i/j{x)dx~^ G r(f2“5), then 

A^s = {fj" -\-vfii) 


and therefore. 


M 3 / 2 ) 


AyS 


f{fi''' + v' fj + V fi') + + vfi) 


dx^/^ 
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On the other side, 


and thus, 


C 


(- 1 / 2 ) 






s = dx ^ 




da:3/2 


= '4’" + V f ^/j' - ^{v f + 


Consequently, 


/: 5 ?f ^ A^s - ^ (^v' f + f + ^/"' + /'^ 


from where we see that the invariance condition fl2^ implies fl24l) . 

Each solution of the projective vector held equation provides a local projective vector held, 
while a global solution of the equation, with the mentioned periodicity condition dehnes a 
global projective vector held. 

It will be proved in next subsection that these projective vector helds can alternatively 
be seen to be the elements generating the stability subalgebra of the point (—l,ndx^) of the 
Virasoro orbit of the corresponding v (see next). 

Let us remark that sometimes v in fl2T|) is replaced by fcn, where A; is a real number,and 
then equation (1251) becomes: 


f + Akvf + 2kv’ f = 0. (25) 

- i.e. / is a solution of the third-order linear diherential equation 

y'" + Akvy'+ 2kv'y = 0. (26) 

3.1 Virasoro algebra and projective vector field equation 

Let Dih_|_(S'^) be the group of orientation preserving diheomorphisms of the circle S^. We 
represent an element of Dih+(S'^) as a diheomorphism <I)(e*^) = where / : M —> 

M is a function such that (a) / G C°°(R), (b) f{x + 27r) = f{x) + 27r, (c) f'{x) > 0. 
Therefore the tangent space TidDih+(S'^) is the set of elements f{x)d/dx G X(S'^) with / a 
periodic diherentiable function such that f'{x) > 0. Not the however that under a change of 
coordinates / does not change as a function but as a vector coordinate. 

The group Diff+(S'^) is endowed with a smooth manifold structure based on the Frechet 
space C°°(S'^). Its Lie algebra g is the real linear space (T(S'^). It is known that the group 
Diff+(S'^) has non-trivial one-dimensional central extensions by (7(1), the Bott-Virasoro group 
Diff_|_(S'^). Recall that a central extension is given by an exact sequence of groups 

1-- (7(1) —A l5iff+(^i) Diff+(^i)-- 1, 
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where j{U{l)) lies in the centre of Diff_|_(S'^). Choosing a normalized (local) section ^ for tt 
we see that o ( 72 ) differs from the product of ^(cxi) and ^{( 72 ) in an element c(ai,a 2 ) of 
j{U{l)). The associativity is equivalent to the cocycle condition 


c(cri o cTa, 0 - 3 ) + c(ai, as) = c(ai, as o as) + c(a 2 , as). 

Changing the section amounts to modify the cocycle c by a coboundary ( (i.e. a cocycle for 
which there exists a map r : Diff+(S'^) —)■ U{1) such that C(c’'i)<^ 2 ) = 'r(ai o as) o (r(ai))“^ o 
(^(^ 2 ))-^). _ 

More explicitly, the Bott-Virasoro group Diff+(S'^) is the central extension of Diff+(S'^) 
by U{1) determined by the Bott cocycle [37] 

c : Diff+(5^) X Diff+(5^) ^ M 


given by 


c(ai,a2) = 


'51 


log(ai o as)' dlog ja^l, 


for ai G Diff+(S'^). This cocycle satisfies the normalization condition c{ai,ai ^) = 0. The 
cocycle can also be written as 


c(ai,a2) = 


'51 


log(a( oas) ciloglaal, 


because using the chain rule we have (ai o as)' = (a( o as) a^ and then, 


1 

2 


'51 


log(ai o as)' dlog ja^l = - log(a( o as) dlog ja^l + - / log(a 2 ) dlog ja^l 


'51 


'Si 


and integration by parts shows that the last term vanishes as a consequence of periodicity. 
Note that 

c(ai o as, as) = ^ / log(ai o as o as)'dlog jag] = ^ / log(a( o as o as) dlog jag] + c(a 2 , as) 


'51 


'SI 


and similarly. 


c(ai,a 2 oa 3 ) = - / log(ai oasoas)'dlog jasoagl = - / log(a( oas oag) dlog jag]+c(ai, as) 


'si 


'si 


from where the cocycle condition follows. If (t, a) denotes the product j(t)^{a) the composi¬ 
tion law in Diff(S'^) is 


(fi, ai) • (fs, CTs) = (ti + t2 + c(ai, as), ai o as). 

The corresponding non-trivial central extension of the Lie algebra X(S'^) by the trivial Lie 
algebra M is called the Virasoro algebra and denoted oir - i.e. we have a central extension 

0--M--oir--X(5i)--0. 
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More specifically, the elements of Dir can be identified with pairs (real number,27r-periodic 
function). In other words, the 27r-periodic function is the component of an element of the 
set which is known to be endowed with a Lie algebra structure. The corresponding 

cocycle ( : X(5'^) x X(5'^) —)■ M is given by 


aXf,X,)= [ fg"dx, 

Js^ 

where Xf = f{x) d/dx and Xg = g{x) d/dx. Obviously, ({Xf,Xg) = —({Xg,Xf) as a simple 
integration by parts shows, because of the periodicity of functions / and g. 

In fact, the corresponding cocycle will be given by 

cP dP 

C{Xf,Xg) = {(^{4>t,^s))\s=t=0 ~ dtds ’ 

where (pt and (ps are the flows of X/ and Xg, respectively. 

Taking into account that 


d 


^(c(0o<^s))|i=o = 9 M{(j)oo'ijj,y){f op,yd\ogp'^ = - I if op,)d\ogpy 


'SI 


si 


and 


d (1 
ds V2 Jsi 


{f oip,)d\ogip'A =- 

/ s=0 ^ 

and similarly for the other term, we hnd that 


'si 


f'dg', 


aXf,Xg) 


\ [ if dg' - g' df) = [ f'g" dx. 
^ Js^ Jsi 


Therefore the commutator in Dir takes the form 


(a,X,),(b,X,) 


(C(W.A',),[W,Ag) 


(y /VXa:, (/s'-s/'jTV 


The dual linear space Dir* can be identihed to the set {(p.,| /a G M,u G 
In fact, an element {g,,vdx^) maps linearly Dir into the set of the real numbers as follows: 


{{fi,vdx‘^),{a,Xf)) = afi+ / f{x)v{x) dx, 

Js^ 

and conversely, each linear map from Dir into M can be represented as such a pair {fi,v dx^). 

The remarkable point is that if ad(a,X/) is the image under the adjoint representation of 
the element {a,Xf) of the Virasoro algebra Dir and ad*^^^^) is the adjoint element, then 


ad(a,x^)(l,'D dx'^) = + 2v f + v' /j . 
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In fact, it follows from the definition 

(ad(a,x^)(/i,'y da;2), (6,Xg)) = {{n,v dx‘^),ad(^a,Xf){b,Xg)) 

= dx^), J'g''dx,[Xf,Xg]^^ 

= h/ /V'+ [ v{fg' - fg)dx = [ (/i/'" -2vf- v'f)g{x)dx, 
Js^ Js^ Js^ 

and then 

ad la,Xf) (h, V dx^) = (0, /"' - 2v f - v' f). 

Hence we see that the stability Lie algebra of the point (—1/2, v dx"^) relative to the action 
dehned by ad* on its dual is given by an Xf G X{S^) such that / is a solution of the third 
order differential equation fl2^ . 

3.2 Projective vector field equation and its structure 

Equation (|2^ is a linear equation and it is therefore invariant under the dilation vector held. 
Lie recipe amounts to introduce a new dependent variable z instead of y in such a way that 
the dilation vector held has the form D = d/dz (i.e. y = Ce^, for any constant C), and then 

y' = Cz'e\ y" = C{z"+ z'^)e\ y'" = C{z"' + dz'z" + z'^)eX (27) 

Note that u = z' = y'/y, no matter of the value of C. The important point is that the 
transformed equation does not depend on z but on its derivatives and the order of the equation 
is reduced by one just introducing the variable u = z' and so we obtain 

u” + ?)uu' + u^ + dvu + 2n' = 0 . (28) 

The purpose of this section is to set up a geometrical formulation of such nonlinear second- 
order ODEs in terms of the projective vector held equation. 

Proposition 1 Let 'ipi o,nd 'ip 2 be two linearly independent solutions of the second order dif¬ 
ferential equation 

fj"+ kv'ilj = d. (29) 

Then, the three-dimensional linear space of solutions of the third order equation 

y'” L-dkvy'+ 2kv'y = d (30) 

is spanned by the functions V’iV’ 2 - 

Proof If "01 and 02 are two linearly independent solutions of the equation fl2^ then taking 
derivatives we obtain 

0j" -\- kvip'i-l- kv' 01 = 0 , -I- fcn 02 + 02 = 0 . 
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Now, if we make use of these two equations, then the following third order derivative 




can be rewritten as follows 

= -{kvi/j'i + kv’'4)i)'ipj — - 'ipi{kv'4)’j + kv''ipj), 

that after simplihcation it becomes 


D^{'4)i'4)j) = -k[2v''ipi'ipj + 4:v{'4)['ipj + ipi'ip'j)]. 


We have therefore obtained 

D^^'ipi'ipj) + Akv{'<p[i)j + V’iV’j-) + 2,kv'{'tpii)j) = 0 , 

what proves that the three functions fij = 'ijji'iljj, i,j = 1, 2, are solutions of fl30|) . Finally, the 
Wronskian of these three functions is given by 

W[^pl, = 2(V’l ^>2 - '^2 'Ip'lf ■ 

Therefore, as is a fundamental set of solutions of the second order equation fl29|) . 

the functions 'ipf, ip 2 and ' 0 iV ’2 are linearly independent and they span the linear space of 
solutions of (130|) . 

□ 

Proposition 2 1. Let f be a solution of the projective vector field equation [2^ . Then the 

function u such that ku = f'/f (where k t)) is a particular solution of the following 
second-order Riccati equation: 

u” + ?)kuu' + k'^u^ + Akv u + 2kv' = 0 . (31) 

2. Suppose that ufix) is a solution of the Riccati equation 

c + kC+ kv = t). (32) 

Then the function u = 2mi is a solution of the second-order Riccati equation ^3l\) . 

Proof: 1. The equation fl26|l is invariant under dilations, and then we can choose an adapted 
coordinate for which dilations generator is {l/k)d/dz (i.e. we use a function 2 ; instead of y 
such that y = e^^), and then using the expressions analogous to fl27)) 

y' = kz' , y" = {kz" + k^R^)e^^ , y"' = {kz"' + Sk^z'z" + kh'^)e^ , (33) 

the differential equation that we obtain does not depend on 2 ; = y'/{ky) but on its derivatives, 
z\ z” and z'", and if we define u = z', it becomes flHT]) . which is a generalization of the nonlinear 
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oscillator equation. Here the coefficients are fixed by the projective vector field equation. In 
other words, the new reduced equation can be rewritten as 

R^{u)+AvR^{u) + 2v' = 0. 

2. First, we note that the derivative of nj + ku\ + /cn = 0 leads to the following second 
order equation 

u'[ + 2kuiu'i + kv' = 0 . (34) 

Now putting u = 2ui in the left hand side of we obtain 

u” + Sfcu u' + k'^u^ + Akv u + 2kv' = 2u'[ + 12 kui u\ + Sk‘^u\ + Skui v + 2 kv\ 

and therefore taking into account fl34p : 

u” + 3kuu' + k'^u^ + Akvu + 2kv' = —2{2kuiu\ + kv') + 12kuiu'^ + 3k‘^u\ + SkuiV + 2k v' 

which symplifying terms reduces to 

u" + 3kuu' + k'^u^ + Akv u + 2kv' = 8/cui (u) + ku\ + kv) 

from where using fl32p on the right hand side we obtain the result of the Proposition. 

This shows a relationship among solutions of ordinary Riccati equations and those of 
related second-order Riccati equations. 

3.3 Global projective vector field and integrable ODEs 

Since the solution space of the projective vector field equation is spanned by 

Span('0^,?/>^,?/>ii/>2), 

an arbitrary solution of the projective vector field equation is given by 

T = H'i/'i + 2i?V'i'02 + Cijjl, (35) 

an arbitrary linear combination of basis vectors. This is periodic when ipi ad '02 are periodic, 
and hence it is a global solution of the equation projective vector field as a consequence of the 
existence of gobal solutions of the Hill equation. This 4/ is called the global projective vector 
field [43] . 

Milne-Pinney equation can be written as a system of two first-order differential equations 
and it turns out to be a Lie system (see e.g. HU)- The same is true for the Hill equation (or 
equivalently for the harmonic oscillator with a time-dependent frequency). As the Vessiot- 
Guldberg Lie algebra of both equations is the same we can determine a mixed superposition 
rule allowing to write the general solution of Pinney equation in terms of two independent 
solutions of Hill equation (see HH): 
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Proposition 3 //V'l ip 2 satisfy Hill’s equation - i.e. they are periodic solutions of 


= 0, i = l,2, (36) 

then the square rooti/j of the global projective vector field, that is, A'lp'l + 2 Btljiip 2 + C'V’i 

is a periodic function satisfying the Milne-Pinney equation mw 

f)” P v{x)f) = (37) 

with a = AC — . 


Note also that we can consider the second order nonlinear differential eqnation [12] 

//''-i/“'+2!.f-| = 0, (38) 

which with the change of variable / = —'ip‘^/2 becomes fl37p . Moreover, taking derivative 
with respect to x at eqnation fl55D we see that a solntion of snch eqnation is a solntion of the 
projective vector £el eqnation fl2^ . 

In a similar way, we can consider the second-order Knmmer-Schwarz eqnation: 


2 / 



-h crp + v(x) 


0 , 


which is a particnlar case of the second-order Gambler eqnation and it was recently analysed 
in |52| from the perspective of Lie theory. It has been proved to be a Lie system associated with 
a Vessiot-Gnldberg Lie algebra isomorphic to sl(2,M), and therefore admitting a nonlinear 
snperposition rnle m. Bnt as in the case of Milne-Pinney eqnation we can also hnd a mixed 
snperposition rnle in terms of solntions of Hill eqnation. More explicitly, the solntion is given 
by 

f(x) = + 2B'ipi'4)2 + (39) 

where ipi and 'ip 2 satisfy the Hill’s eqnation and B^ = AC — aW~^, where W is the Wronskian 
determinant of both solntions, W = hP['^i,'^ 2 ] [52]. 


It has also been proved in [63] that if fji and V '2 satisfy Hill’s eqnation, with '02(a^o) 7^ 0, 
and W denotes the Wronskian, W = W[fji,'ip 2 \i then for each nonzero real nnmber m, snch 
that 0 7 ^ m 7 ^ 1, the fnnction 





c 



l/m 


satisfies the differential eqnation 


y” + v{x)y = c 


(^iV’s 


\m—2 


y 


2 m—1 
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There is a clear difference with respect to the preceding case of Milne-Pinney equation because 
now the right-hand side of equation depends on the functions ipi and ip 2 - 

This result was generalised in [M] where it is shown that if 'i/ji and "02 are functionally 
independent solutions of the linear homogeneous differential equation 

y” + r{x) y' + q{x)y = 0, 


then the function 

is a solution of the Reid type equation 

y" + r{x) y' + q{x) y = AB{m - l)(V’iV’ 2 )”' 

where W is as before the Wronskian of the two functions. 

Moreover, if we set "0 = (' 0 i' 02 )^'^^ then -0 satisfies Thomas equation [ 68 ] 

7/2 1 

y"+ r(x)y'+ kq(x)y= (1-m - kW^y^-^\ kl = 1, (40) 

y 4 

where W denotes the Wronskian of the two solutions. 

We can recover reduced Gambier equation (or PainlevGGambier XXVII equation) [H] for 
special values of m. 

4 Invariants, prolongations and Riccati equation 

For a geometric approach to a partial differential equation in m independent variables and 
one dependent variable v, we must consider the space M x V = {(x, u) | x G M,u G V}, 
where M = R™' and V = R. Suppose that G is a Lie group acting on some open subset 
N O M X V. Then the transformation by 5 ^ G G is 

fi'(x,M) = ( x , m ), geG, 

and a hypersurface u = 'u(x) in N is transformed into another one, u = 'u(x). 

An infinitesimal transformation is given by 

X = X -I- ef(x, m) -I- 'd(e^), 

u = u + e g{x,u) + 

wich can be understood as the flow of the vector field in M x V 

We restrict ourselves to the case of ordinary differential equations (i.e. m = 1) and then 
M = R. A (may be local) smooth section a for the projection tt : M x V ^ M defines a 
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smooth function u = u{x), by means of (t{x) = {x,u{x)) and then for each natural number k 
induces a function = pr’^^^'u(a;), called the k-th. prolongation of u, where 

pr^^^M : R —^ R^+^ 


is the curve whose components are the derivatives of u of orders from 0 to k. The total 
space M x C R*^+^, the coordinates of which represent the independent variable x, the 
dependent variable u and the derivatives of u to order n is called the n-th order jet space of 
the underlying space M x V, sometimes denoted 

Similarly, a vector held m N <Z M xV given in local coordinates by 

X = f{x,u)— +g{x,u)—, (42) 

admits a fc-order prolongation. We are only interested in the hrst-order prolongation [58]. 
Note that if we consider the corresponding inhnitesimal transformation flTT|) . then 

du 

{^9 UxDf) 9x {9u fx'jUx fuUx-i (43) 

where Qx = dg/dx and g^ = dg/du, and similarly for and fu, with D being given by 


^ d d 

D — — -h 

ox ou 


It follows directly from 

M _ d{u + eg + ^(e^)) _ Ux + {gx + guUx)e + 'd(e^) 
dx d{x + ef + d{e^)) 1 + (/^ +/«Mx)e +'^(e^) 

= Ux + e {Dg - UxDf) = Ux + [gx + {Qu - fx)ux - fuul]e + '^(e^). 


This provides us with the following dehnition of hrst-order prolongation of X|58]: 

Definition 3 The prolongation pr(^)(X) of the vector field X G X(iV) is the vector field in 
given by 

prW(X) = = f{x,u)— + g{x,u)— + {Dg - UxDf)—. (44) 

These prolongations will play a relevant role in the search for symmetries of diherential equa¬ 
tions as it is shown in next subsection where the hrst-order Riccati equation is used as an 
example. 
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4.1 Infinitesimal symmetries of standard Riccati equation and pro¬ 
jective vector field equation 

The first-order differential equation for the function u 

Ux = V? + v{x), (45) 

where v is now a given function, usually called standard Riccati equation, is a particular case 
of the general Riccati equation: 

Ux = a 2 {x)v? + ai{x)u + aQ{x) ^ (46) 

for 02 = 1, Oi = 0 and a^^x) = v{x). 

The solutions of such non-autonomous differential equation (145)) are the integral curves of 
the vector held in given by 

A-= I-+ + .„(.))!-, (47) 

An inhnitesimal Lie symmetry of such differential equation is represented by a vector held 

such that there exists a function r satisfying 

[Y,X]=tX, (48) 

because then the solutions of the diherential equation are transformed under the how of Y 
into solutions, up to a reparametrization. 

Taking into account that 

[fdx + gdu, dx + {u^ + v)du] = -{Xf)dx + {fvx + 2 gu - dxQ - (u^ v)dug) 

we hnd that the symmetry condition implies that r = —fx{x), and that / and h are related 
as follows: 

8 x 9 + + v)dug - fvx - 2gu = {u^ + v) fx . (49) 

We can alternatively consider the given diherential equation (1451) as dehning a 2-dimensional 
submanifold in R^, with local coordinates {x,v,w), dehned by zero level set of the constant 
rank map 0 : R^ ^ R, (j){x, u,w) = w — — v{x), that is, 

S = 0“^(O) = {{x,u,w) \ w — — v{x) = 0}. (50) 

The hrst-order prolongation of Y is given by 

Y^8 = f(^x)—+g{x,u)— + {gx{x,u) -wfx{x))— (51) 
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and when we consider the symmetry condition, which is but a tangency condition, = 0 , 

we hnd 

+ v{x)))i^=u 2 +y(x) = 0 , (52) 

and more explicitly, 

d^g{x, u) + {wdug{x, u) - /x(a:)w)|^=„ 2 +^(^) = f{x)v^{x) + 2g{x, u) u , (53) 

which reproduces (09]). 

We can now establish the following relationship among symmetries of the standard Riccati 
equation and solutions of the projective vector field equation. 

Proposition 4 The standard Riccati equation ( 03] ) remains invariant with respect to the first 
prolongation pr^^^F of the vector field 

Y = nx)^ + g{x,u)^^ 

provided that f satisfies the projective vector field equation. 

Proof. Once that / is fixed in f09|) we can try to determine an affine function g (i.e. of 
the form g{x,u) = k{x)u + h{x)), satisfying such symmetry condition. Then the functions k 
and h must satisfy: 

k^{x)u + h^{x) + (u^ + n(a;))(fc(a:) - ^(a;)) = f{x)v^{x) + 2u{k{x)u + h{x)), 
and therefore 

u^{k{x) + fx{x)) + u{2h{x) - k^{x)) + f{x)v^{x) + v{x){f^{x) - k{x)) - h^{x) = 0. 

This shows hrst that we must choose k{x) = —fx{x) - i.e. the function g is such that 

g{x,u) = -fxix)u + h{x), 

and using such an expression for g{x,u) in the symmetry condition f09|) we hnd 

u{fxx{x) + 2h) + f{x) Vx{x) + 2v{x)fx{x) - hx{x) = 0 , 

which shows that we should choose h such that h = —^fxx{x), and then replacing h by this 
value in the preceding equation we hnd that / must be such that 

^fxxx{x) + 2v{x)f,:{x) + f Vg;{x) = 0 , 

and hence as indicated in |13] / is a solution of the projective vector held equation (1241) . 

Remark: One would obtain the same result if one starts from the diherential equation 
u' = a 2 {x)u^ + ai{x)u + ao{x). In this case u must be expressed in terms of 02 , Oi and Oq and 
their derivatives. 
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4.2 Integrals of motion and other dynamical features 

A constant of the motion or symply an integral of the motion for a system of ODE’s 


dVi 

dx 


■ ■ ,yn) i 


1 , • • • ,n 


( 54 ) 


is a non-constant differentiable fnnction <h(x, yi, • • • , ?/„) that retains a constant value on any 
integral curve of the system. This means its derivative with respect to x vanishes on the 
solution curves: 


dx 


^ 9 $ dyt 

dx dyi dx 


£>[4] = 0, 


(55) 


where 


D 


A 

dx 


i 


A 

dyi 


is called the material or total derivative. For an autonomous system this reduces to 


i 



(56) 


where 


D = Y,^‘ 

i 


A 

dyi 


is just the vector held associated with the given autonomous system. The a:-independent 
integrals of motion are usually called hrst-integrals. 

Higher order differential equations can be written as an associated system of hrst-order 
differential equations and then the constants of motion for such system are called constants 
of motion for the higher-order differential equation. As an instance the projective vector held 
equation fl2T)) can be written as the system 


( dy 



< — 
dx 
da 

< dx 


w 

a 

—ivw — 2v'y 


and then the function ^{x,w,a) = 2v{x)y‘^ — + ya is & constant of the motion, because 

the vector held 


is such that X<h 


X 


A A 

dx dy 


d 

a— - (Avw -|- 2v'y) 

dw 


da 


0. The constant of motion can also be rewritten as 


® ( 2 :, y, y\ y") = 2v{x)y^ - -y'^ + yy 


n 


22 




which is a constant of the motion for the projective vector field eqnation (Ell). This shows 
that a nowhere vanishing fnnction is solntion of the projective vector field eqnation fl2l)) if 
and only if there exist a constant C snch that the fnnction is solntion of the second-order 
differential eqnation 


y 


n 


2v{x)y 


^'2 

y 


With an analogons procednre one can compnte the first-integrals of similar type for a 
third-order differential eqnation. For example, the first-integral of the stationary Calogero- 
Degasperis-Ibragimov-Shabat eqnation [25] . 


y"" + 3y^y" + 9yy'^ + 3y%' = 0 


(57) 


has associated a vector field 


d d 

X = w— -I- a— - {3y‘^a + 9yw‘^ -|- 3y^w) 

ay ow 


A 

da 


and the fnnction 


<F(|/, w, a) = 2ay + 6y^w + y^ — 


is snch that = 0 and therefore is constant of motion. The corresponding fnnction given 
by 

y', v") = y{2y'' + (6yy' + y^y) - y 

/(2/.. + (fifU + f‘)f) - fJ, 

is a particular case of the previous one with v = {3/2)yy' + (1/4)?/^. The remarkable point 
here is that one can interpret equation fl57|) in terms of the stabilizer set of the coadjoint 
action. The function y = f{x) is a solution of (l57|) if and only if 


ad W (3//,r + (1/2) f'^)dx‘^ = 0. 

dx 


5 Second-order Riccati equation, Painleve II and higher 
Painleve type equations 


Consider once again the linear third-order equation fl2l|) and the corresponding second-order 
Riccati equation (12811 . 

Let us define for each function / the functions 



f 



where w is an arbitrary but fixed function. Then, 


u 


/ 

1 


d 

dx 



r 

f 



U2 - -Ui - W. 


(58) 
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(59) 


Moreover, 


rff/ £11 £! 

f J J J f , / 

^2 = -^- p -^ 

and if we use the second equation in fl58|) to write 

f" 1 2 

— = U2 + -U^-W, 

we see that when / satishes equation fl 2 ^ and therefore 

f" 


f 


= —4vui — 2v', 


then 


u '2 = —4v Ui — 2v' — Ml ^M 2 + — uiu[ + w' 

u '2 = —4mM l — 2v' + Ui{w — M 2 ) — Ml ^m'^ + + w'. 


that is, 

Consequently, such functions Mi and M 2 satisfy the system of differential equations 


m; 


u'2 


1 2 

U2 - -ui 


w 


= 2{w — 2v)ui — 2 miM 2 — 2v' + w'. 


(60) 


Conversely, if mi and M 2 are solution of this system then f{x) = exp mi(C) dQ is ^ solution 
of (Ell), because, by dehnition of /, 


11 . 

f 


Ml, 


r 

/ 


d 

dx 



+ 



M]^ + 


U2 + -Ml - W, 


and 


r 

/ 


d 

dx 



+ 


rr 

p 


u '2 + Uiu'^ — w' + Ui ^M 2 + -u\ — = u '2 — w' + 2 Mi(m 2 — m;). 


from here we see that 


/"' + 4m /' + 2m' / 

7 


= u '2 — w' + 2 Mi(m 2 


ta) + 4 mmi + 2 m' 


and the second condition ( 16 ( 11 ) shows that / is solution of (l 2 ll) 

The very essence of the second-order Riccati equation fl28|l is that it is obtained as a 
system of two coupled differential equations. One starts from a hrst-order Riccati equation 
u'l = M 2 — \u\ —w, where w is an arbitrary differentiable function of the independent variable 
X and then couple Mi to M 2 through a linear differential equation in M 2 involving mi. 
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We wish to express these two equations in terms of U 2 - It is clear from the second equation 


that 


Ml = 


1 Mn + 2v' — w' 


2 w — 2v — U 2 

and if we define a new function K{x) hy K = ±(m 2 — w + 2v) we can write 

1 K' 

“ ~2'K' 


and if we reexpress the hrst equation of fl60l) in terms of K and substitute U 2 hy w — 2u±K 
we obtain 


K" = -- 


4 K 


+ AvK T 2K\ 


(61) 


We can put it into a more amenable form if we substitute K = —1/p, and now we obtain 


p" = - - Avp =F 2. 

4 p 


/2 


(62) 


This is a reduced Gambier equation or G5 equation in Gambier’s classihcation m 

Remark (a) If we assume n to be a constant such that n = ^, then the second-order 


(63) 


Riccati equation fl25]) boils down to the modihed Emden equation 

u” + ?>uu' + + au = 0, 

which therefore can be obtained from the pair of differential equations 


u'l = U2 — \u{ — W, 

u '2 = {2w — a)ui — 2 uiU 2 + w'. 


(b) Suppose we take u = 2\k'/\I'. Then, the second-order (projective) Riccati equation 
can be related to the family of Milne-Pinney equations 


T" + vA/ = 


a 


(64) 


where a is some constant. In fact, such a family is described by the third order differential 
equation 

(T^T^ + mT^)' = 0 . 

This is a differential equation invariant under dilations and Lie recipe for reduction 
amounts to define a new variable such that 4/ = (i.e. u = C = 24/'/4/), and then 

replacing 

T' 1 T" 1 , 1 , T'" 1 „ 3 , 1 o 

~r ~ ~r ~ 0 “ + 7 "^ > ~ 7 "^ + 7 "^“ + -I 




in the equation 




-h 3T^4>'4>" -F 4 mT^4>' -F = 0, 




r2iTr^ iTr^^ 


3iTf' I 
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we obtain the second-order differential eqnation fj28l) . This immediately yields onr resnlt. □ 


The projective vector field eqnation fl2Tl) admits a Lax formnlation in the following sense. 
Given a fnnction / define a pair of matrices P and Q as follows: 


ff/2 -J^viOfiOdC 
-f -f/2 


Q = 


0 V 
-1 0 


for which 


QP = 


-vf 


-vf'/2 


PO- (T'^^OfiOdC vf/2 

^ V /V2 -vf 


-f/2 r^ofiodcj^ 

Therefore, the matrices P and Q are a Lax pair, that is, 

P'+[Q,P] = 0, 

if and only if / is a solntion of the differential eqnation 

f + 2vf + 2 rv{Of{OdC = 0, 


and taking derivatives we find the projective vector field eqnation 

The corresponding Lax formnlation for the projective second-order Riccati eqnation 
will be given by the same matrix Q and 

< J/-I 


(65) 


_ I ^ued “(f) v(()u(()ed “(d)‘^d(i(j’ 


P= 2 




l.yeR“(C)<iC 


and then the matrices P and Q are a Lax pair if and only if n is a solntion of the integro- 
differential eqnation 


+ + + 2 / n(C)M(C)e^ dC = 0. 


Now, taking derivative with respect to x we find 

(u" + 2uu' + + u(u' + + 2n)e^" = 0, 

and simplifying the common factor we obtain the second-order Riccati eqnation 0281) . 

Remark The eqnation fl65l) has a nice geometric interpretation in terms of the coadjoint 
action of loop algebra [62l |65]. The loop gronp C'°°(S'^, G) is the gronp of smooth fnnctions 
on the circle with values on real semi-simple Lie group G. Its Lie algebra G°°{S^,g) has 
a central extension provided by the Kac-Moody cocycle uj{A,B) = < A{x), > dx. 

The coadjoint representation of the extended loop algebra is given by 

ad*{P,a){Q,h)= ([P,g] + 6^^,0), 

for any P,Qe G°°{S^,g) and a, 6 G M. The equation fl6^ is the stabilizer set of the coadjoint 
orbit confined to hyperplane b = —1. Thus once again we obtain equation which can be 
interpreted in terms of stabilizer orbit. 
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5.1 Painleve II, Backlund transformation and second-order Riccati 
equation 

We wish to discuss integrable class of Painleve II equation; in other words, we are interested 
in the rational solutions for integer valued parameter a of the Painleve II equation and their 
explicit characterization in terms of the Airy function. It is worth to note that half-integer 
valued parameter of the Riccati equation are also characterized by the Airy function. Here 
the Painleve II equation (PII) is an ordinary nonlinear second-order differential equation with 
a parameter a, 

u” = 2u^ + XU + a. (66) 

This equation has exactly one rational solution for a being an arbitrary integer and has no 
rational solution if a is not an integer [T8l [28l [29] . It admits a Backlund transformation 
(u, a) I—)■ (n, —a), which clearly maps rational solutions into rational solutions. 

Consider now the Airy differential equation 

^l;” + x^p = 0. (67) 

It is clear that the Airy differential equation is a particular instance of the Hill’s equation 
fl36|l for the special choice v{x) = x. 

Proposition 5 1.- If'ip is a solution of the Airy equation 

fj” + xip = t), 

then the function ui = fj'/ip satisfies the Riccati equation u' + + x = t). 

2.- The corresponding second-order Riccati equation for the function u 
Painleve II equation 

2 

u” = -f- 2a; M — A, 

which after a rescaling becomes the Painleve II equation 


= A Ml is the 
( 68 ) 


Proof. 1.- It suffices to use the expression (IT^ for reduction recipe with n = 2, ao{x) = x 
and ai(a;) = 1 (i.e. R^{u) xR^{u) = 0), and expressions (ITTH . then ui = ip’/ip is a solution 
of 

m' -I- -I- a; = 0. 

2 

2. Note that the function u = Xui satishes the differential equation u' = — A a; and 

therefore, deriving with respect to x we hnd that the second-order Riccati differential equation 
satished by the function m = A Wi is 


u" = — 


2u 


u 

T 


A X 


- 2m3 

A = —- + 2xu 
A^ 


A. 


If we now consider a new independent variable x = fix with fi ^ 0, the differential equation 
is transformed into 


d?u 2u^ 


X A 



and by choosing the parameters A and /i as /i^ = 2 and A = /i ^ we see that the differential 
equation becomes a Painleve II type differential equation: 


d?u 

dx‘^ 


2u^ + XU — 


1 

2 ‘ 


This proves that the second order Riccati equation can be brought back to ‘standard Painleve 
IP by an appropriate scaling. □ 

Under the map x —)■ —x, it takes the form of a similar equation but with the opposite 
sign, that is, i.e. if is a solution of the Airy differential equation 


'ijj" — = 0, 


then Ml = AV’VV’ satisfies the differential equation 

2 

u” = —u^ — 2xu — X (69) 

A"* 

that can be transformed into a Painleve II equation with an appropriate change of independent 
variable as before. 

Remark The Painleve transcendents (P-II - P-VI) possess Bdcklund transformations which 
map solutions of a given Painleve equation into solutions of the same Painleve equation, but 
with different values of the parameters. Therefore two Painleve equations for a = 2 and 
a = —2 are connected by Backlund transformations. 

Remark Considering ip"±{x/2)ip = 0 we can also relate Painleve II equations for parameters 
a = =f 2 and these are connected by the Backlund transformations. 

Remark li y = f{x) is a solution of the differential equation fl2T)) . then for a given real 
number /3 the function y = f{x) is a solution of the differential equation 

y"' + Avy' + 2v' y + Qfdy y' = (70) 


with V = V — /3 f. Now if the function / is positive and define the positive function if by 
f = then as 

f = 2'ip'ilj', f" = 2-^'^ + 2'ip fj", f" = if” + 2-^ tp"', 


we see that p) is then a solution of 


2y y'” + Qy' y” + %vyy' + 2v' y^ + 12/3 y^y' 


0 , 


which can be rewritten as 

{y"y^+ vy^ +^y^y = 0, 

i.e. there exists a constant a such that is a solution of the differential equation 

y” + vy + Py^ = 

y6 

For the particular case v = x/2, that means that pj is & solution of the Ermakov-Painleve II 
equation. Conte showed how one can transform it to Painleve II (see e.g. [66] ) 
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Remark Let us briefly describe the connection between Painleve II hierarchy and our ap¬ 
proach. It is readily clear that the projective vector held equation is the stability algebra 
of Virasoro orbit. In other words, instead of considering the symplectic structure dehned by 
Oi = dx and the Hamiltonian function i.e. 





A 

dx 



UUx + u 


XXX I 


we consider the second Hamiltonian structure of the KdV equation 


O2 — 0^. + -udx + ( 71 ) 

together with the Hamiltonian density Using “frozen Lie-Poisson structure” we can 

dehne the hrst Hamiltonian structure of the KdV equation too. This satishes famous Lenard 
scheme 


OxHu+i = 


2 1 

dx + —udx + -^Ux 1 'Hr, 


(72) 


where 


Hi = u, 




2 

X'f 


are the conserved densities. The mKdV hierarchy is obtained from the KdV hierarchy through 
Miura map u = Vx — v'^. The second Painleve hierarchy is given recursively by N. Joshi from 
the modified KdV hierarchy |16] 


P"; («, = ( T + ) Jn{u,r - U^) - X U - fin = 0, 


(73) 


where /3n are constants and J7„ is the operator defined by the first and second Hamiltonian 
structures of the KdV equation 


dx>Jn+i{u) — {d^ + 4:udx -l- 2ux)J^n{u) 


with Ji{u) = u. 

Lemma 1 Let'01 und 02 two linearly independent solutions of Hill’s equation. Then, 

(a) The differential equations 

+ IQvy” + lOuV + (Qu^ + ?>v'')y = 0 (74) 

trace out a four-dimensional space of solutions spanned by {0f, 0 f 02 , f’i}- 

(b) The differential equations 

y^” 20vy"' 30v'y'' -|- 18v"y' -f- QAv'^y' -|- -|- -\-64:Vv'y = 0 (75) 

traces out a five-dimensional space of solutions spanned by {0f, ' 0 f' 02 , '^ 2 }- 
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Proof. By a direct lengthy computation. □ 

Using the standard projective transformation u = y'/y we obtain the third-order and 
fourth-order Riccati equations associated to these equations, respectively, which are, according 
to flT^ given by 

u"' -|- 4m u" -|- u' + lOn u' + u'^ + lOn + lOn' u -|- -|- 3v" = 0 (76) 

and 

-I- 5m v'" -|- IOm'm" -|- 15mm'^ -|- IOm^m" -|- IOm^m' + + 20m(m" -|- 3mm' -|- m^) 

-I- 30m'(m' -|- m^) -|- 18m"m -|- 64m^m -|- 4m'" -|- 64mm' = 0. 

Note that if we put m = 0 then, all these equations form a Riccati hierarchy. 

5.2 Higher Riccati equations and higher-order Painleve class sys¬ 
tems 

One can easily check that the third-order Riccati equation fl76|) can be transformed by setting 
M = M^ into a special case of the Chazy equation XII, given by 

m'" + 10mm" + 9m'2 + 36m^m' + 20m^ = 0. (77) 

Recently Ablowitz et al DEI studied a general class of Chazy equation, dehned as 

m'" — 2mm" -|- 3m'^ = q;(6m' — m^)^. (78) 

where a is a real number. The particular case a = 1/16 was mentioned in [2T] and that of 

4 

^ ~ 36-n^ ’ 

m'" — 2mm" -|- 3m'^ =-5'(6m^ — u^y. (79) 

36 — 

has a singled value general solution. More explicitly, this equation was hrst written down and 
solved by Chazy dsidn] and is known today as the generalized Chazy equation. Clarkson 
and Olver showed that a necessary condition for the equation (52) to possess the Painleve 
property is that the coefficient must be a = with 1 < n G N, provided that n y Q. 

It has been further shown in Cl that the cases n = 2,3,4 and 5,correspond to the dihedral 
triangle, tetrahedral, octahedral and icosahedral symmetry classes. 

Proposition 6 If u is a solution of the third Riccati equation 

m'" + 3m'2 + 4mm" + bM^M' + M^ = 0, (80) 

then, using x = —x as independent variable, the function m = 2 mi satisfies 

m'"-2mm" + 3m' 2 = i(6M'-M2)2. (81) 

8 
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Proof. In fact, if n is a solution flSOll then u = 2 m is solution of 

u'” + + 2uu' + -u^ u' + = 0, 

2 2 8 

and with the mentioned change of independent variable, 

u”' — -v'^ — 2uu' + -u^ u' — = 0, 

2 2 8 

which can be rewritten as in fISTD . □ 

Remark The Chazy IV equation 

u'” = -3mm" - 3m'2 - ?>u^u' (82) 

is a derivative of second-order Riccati equation or second member of the Riccati chain gven 
in ([II]) with k = l. 

We are able to construct fourth-order equations of the Painleve class family, derived by 
Bureau [5]. The Painleve classihcation of the class of differential equations of the 4-th order 
hrst and second degree was studied by Cosgrove [20] [2l]. The subcase which will be relevant 
here is the Bureau symbol PI. 

Cosgrove presented the results of the Painleve classihcation for fourth-order differential 
equations where the Bureau symbol is PI. He gave a long list of the equations F-VII - F- 
XVIII in this category. Six equations, denoted by F-I, F-II, ..., F-VI, have Bureau symbol 
P2. It is worth to note that all the cases with symbols P3 and P4 were found to violate a 
standard Painleve test, they admit non-integer resonances. 

We derive the following different equations from the list of Cosgrove on the Painleve 
classihcation of fourth order equations with Bureau symbol PI. 


Proposition 7 The following two equations follow from the higher Riccati equations 


F-XII 

F-XVI 


m(™ = —4mm"' — 6m^m" — 4m^m' — 12mm'^ — IOm'm" 

^{iv ^ -^uu”' - IOm'mm" - 15mm'2 - IOm^m" - IOm^m' - 

+ A{x){u"' -I- 4mm" -|- 3m'^ -|- 6m^m' -|- m'^) -|- B{x){u'' -|- 3mm' -|- u^{x)) 
+ clx)^ + u') + D{x)u + E{x) = 0. 


Proof A) The F-XII fourth order equations with Bureau symbol PI follows directly 
from the expressions flT^ for the fourth order Riccati R‘^{u) = 0 and the third-order Riccati 
R^{u) = 0 equations (also known as Burgers higher-order hows) by means of the relation 


F-XII : R\u)-uR%u) = 0. 


B) The F-XVI fourth order equations with Bureau symbol PI is the combination of all 
higher order Riccati equations. □ 
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5.3 Second degree Painleve II equation 

We have seen in the earlier section how Painleve II is connected to second-order Riccati or 
projective vector held equation. In this section we carry out this investigation further to 
incorporate the second degree Painleve II equation. 

The Hamiltonian of the standard Painleve II equation u” = 2u^ x m -|- a is given by 

H{x, u,w) = Y - ^ - (^« + 0 (83) 

where u and w stand for coordinate and momentum. The Hamiltonian equations of motion 
yield a set of Riccati equations 


u = 


w = 


dH „ X 

= W-U^--, 

2 1 


(84) 


^ =2uw + a + ^ 
ou 2 


This system was studied by Morales [M] who proved that for a G Z the system is not 
integrable by means of rational hrst integrals. 

Taking derivative with respect to x at the hrst equation, using for w the value obtained 
form it, i.e. w = u' + v? + and replacing w' by the value given by the second equation we 
obtain 

u” = 2u {u' + + - — 2 mm ' — - = 2v? + ux -\- a, 


and therefore u satishes the Painleve H type diherential equation 

Now, if u{x) and w{x) are solutions of the system flHT)) . the function h{x) dehned by 
h{x) = H{x,u{x),v{x)) is such that 


dH 


w 


w 


h'{x) = —{x,u{x),v{x)) = h"{x) = -y = -^ {2uw + a + 

and then one easily check that 


(h"(x))^ -I- 4 {h'{x))^ + 2 h'{x) {x h'{x) — h{x)) — - I a + 


= 0 . 

Conversely, if a function h{x) satishes the relation flHSD . then the functions 

1 


1 h"{x) f 1 


4h'{x) ’ 


w{x) = —2 h'{x) 


are solutions of the system 

Observe that the right hand side of the expression 


(85) 


X 


2h' = w = u' + + —. 


( 86 ) 
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is the term appearing from the reduction of the second-order linear 
y"+^y = 0 when considering Lie recipe of invariance under dilations, i.e. 
Hill’s equation 


y” + (2 h'{x) + I j 1 / = 0, 


differential equation 
u = y'/y. Therefore, 

(87) 


leads by reduction from dilation symmetry to the preceding equation. 

Note that the solution u can be found directly from the function h 
the second expression in fjSTll : 


u = 


2h" + a + I 
4/d 


the expression using 


( 88 ) 


Thus we can extract several important information about the Painleve II equation from 
the second-order projective Riccati of projective vector field equation. 


6 Outlook 

The Riccati equation plays a very important role in mathematical physics and dynamical sys¬ 
tems. The second and higher-order Riccati equations also play a big role in integrable ODEs, 
Painleve and Chazy equations. It has been shown [8] that the use of geometrical techniques 
to deal with the elements of the Riccati equation is very efficient to unveil some previously 
hidden aspects of such equations. The novelty of this article is to study (higher) Riccati equa¬ 
tions and various other connected integrable ODEs using coadjoint orbit method of Virasoro 
algebra. In particular, we have given its geometric description using the stabilizer set of the 
Virasoro orbit or projective vector held equation. It would be interesting to extend our study 
to coupled Riccati equations, how they are connected to stabilizer set of the extended Vira¬ 
soro orbit or superconformal orbit. The application of these differential geometric methods 
to deal with multicomponent systems and their integrability is a very interesting subject to 
be studied. Hence there are several interesting issues connected to this paper ought to be 
addressed in future, we have only hit the tip of the iceberg. We hope to answer some of these 
questions in forthcoming papers. 
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